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ABSTRACT 



A theoretical and experimental investigation of the stability of 
developing laminar flow in a circular pipe subjected to small, rota- 
tionally symmetric (torsional) disturbances is presented. Analytically, 
developing mean velocity profiles were represented in polynomial series 
form and subsequent formulations yielded determinations on stability 
which reduce to the classical conclusions for the case of fully developed 
flow. 

Experimentally, torsional disturbances were imposed upon a developing 
laminar pipe flow by sinusoidally rotating a 1*5 inch (1 diameter) length 
of pipe in the entrance region of a carefully constructed wind tunnel 
facility. For flows with maximum Reynolds numbers of 8500, maximum 
amplitude and frequency of oscillation were % inch and 30 Hz., respec- 
tively. Hot wire anemometer measurements verified the stability of the 
developing flow field. Moreover, measured wave speeds and decay factors 
compare favorably to calculated values. 
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NOMENCLATURE 



a^, a^, a^, Constants, coefficients of the power series 

solution of the disturbance amplitude function 

^0* ^1’ ^2^ ^3 Constants, coefficients of the polynomial 

approximation of the velocity profile 



C 

r 


Dimensionless wave speed 


C. 

1 


Amplification factor 


d 


Inside diameter of the pipe 


P 


Pressure 


r 


Radius of the pipe 


R 


Reynolds number (R-^) 


U 


Mean velocity in the pipe 


U 

max 


Maximum velocity in the pipe 


U 

o 


Dimensionless velocity profile in the entry 
region, referred to U 


u 


Axial component of velocity 


w 


Tangential component of velocity 


z 


Axial coordinate 


oc 


Disturbance wave number 


l) 


Kinematic viscosity 


y 


Stream function 
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I. INTRODUCTION 



The purpose of the work reported here was the study of the evolution 
of small, rotational ly symmetric disturbances introduced into the laminar 
transition region (the hydrodynamic entrance length) of Poiseuille pipe 
flow. 

Since the pioneering work in the field of hydrodynamic stability 
conducted by 0. Reynolds in the 1880' s, the stability of fully developed 
pipe flow has been the subject of many theoretical and experimental 
studies. All of the analyses up to the present concerning the stability 
of fully developed pipe flow show that such a flow is stable to small 
disturbances • 

A review of the theoretical and experimental determinations con- 
cerning the transition phenomenon leads one to the following considera- 
tions: (1) transition impends because the disturbances introduced 

contain some initially large disturbances so that small disturbances 
are not significant regardless of when or where they are introduced; 
or, (2) initially small disturbances, either axially-symme trie , rota- 
tionally symmetric or three-dimensional, grow in size along a favorable 
length of conduit, namely the entrance length, via convective inter- 
action while the boundary layer is rapidly developing. Consequently, 
these small disturbances, which would not have grown in a fully developed 
region, now acquire a magnitude which eventually renders the fully de- 
veloped flow unstable. 

A theoretical analysis in support of the above considerations has 
been presented by Tatsumi [1]. His conclusions demonstrated that small 
disturbances magnified in the entrance region of a pipe flow field above 
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a minimum Reynolds number of 195400^ do cause transition after a passage 
of 17 tube radiic The significance of his study is not the determination 
of a value for the critical Reynolds number^ but that there is a pre- 
diction of ins tability <, 

Presently, the case of developing axially-symmetric flow is treated. 
Results are shown to be reducible to the determinations of Pekeris [5] 
for fully developed pipe flow. Additionally, experimental verification 
of calculated wave speeds and amplification factors is presented. 
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II. THEORETICAL ANALYSIS 



A. CYLINDRICAL ORR-SOMMERFELD EQUATIONS 

To describe a general disturbance in mathematical form requires the 
formulation of a three-dimensional (r, 0, z) vector function of time, t. 
The generally accepted treatment for small disturbances is to view them 
as a Fourier series of a large number of elementary disturbances, each 
of which can be represented in time and space by a finite number of para 
meters. Indeed, this is the foundation of the small disturbance theory 
that the disturbance is a linear combination of spatially periodic 
components which grow or decay with time. Such a component is character 
ized by the stream function, 

F(r) Exp 

By further contending that the imposed disturbance be rotationally sym- 
metric, then the 0 dependence is eliminated, and the representation 
reduces to, 

Exp[6oc(Ct-^^ (2) 

Therein OC is a real number - the wave number of the disturbance in the 
fluid. C is a complex number, of the imaginary part 

of C, C^, determines whether amplification (C^ <C 0) or decay 3^ 0) 
occurs, whereas is the wave speed of the disturbance. Experimentally 
such a disturbance is difficult to generate, since, as formulated, it 
should be imposed upon the entire length of the conduit over the entire 
cross-section of the conduit simultaneously, and very periodically in z 
as well. 
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Analytically, the disturbance described in (2) is thought of as 
being imposed upon a steady mean flow characterized by a stream func- 



tion, 



Vi = r 



The equations of motion for incompressible flow symmetric about the 
axis of a pipe are: 






_L AT ^ U.^ — — L -f- 2^ 



■ji* '''5^ ' 



(4b) 



M: 



AT - c) 



(4c) 



\7^r 






■-1. i' 



h “y?- 



(4d) 



It is convenient to introduce the non-dimensional coordinates. 



z 



j ? = 



-to o'- P 
5 ^ - V D P ^ ? 



d. J ^ - <i, J ^ - d J ' ^ p 

and define R, the Reynolds number as R = With these substituted, 

the system of equations (4) becomes: 



, at' b/u' , 








X- a'lu'- _ 1.1^ 9'^a' 



(5b) 
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4 - 



^ +- 4^ = o 






(5c) 



Defining U = and V = - ^ , it is noted that with the dependence 

on 0 dropped, that there are two possible types of symmetrical distur- 
bances which may be generated and which are not mathematically coupled. 
First, there exist the torsional disturbances, which are concentric 
with respect to the pipe axis; secondly, there exist the meridional 
disturbances, which occur in planes passing through, and parallel to, 
the axis of the pipe. Treatment of each case is necessary due to the 
aforementioned difficulty in experimentally imposing the mathematically 
exact disturbance upon a flow regime. It is assumed that. 



= f (r) Exp[c'^(Ct -2^ 



(6a) 



SI - <5(0 EXp[(:o( (ct-2-^ 



(6b) 



These relations, together with the definitions of velocity components 
in terms of the stream functions, lead to the following uncoupled sys- 
tem of exact equations of motion: 



+ - If IpO’-rJ + 



(7a) 



(7b) 



where 






— 1 ^ 

h ^y- 



-f 



As is composed of the main flow with the disturbance superimposed 

upon it, it can be described as. 
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(8) 



T = % t r' 

where refers to the main flow and ^ is as defined in (6a), 

Substituting (8) into the system (7), dropping squares and products 
of and , as each is assumed small p leads to the following 

uncoupled system of linearized equations: 






p" 








kJL 



+■ ^ 



r 

L^r 



^JL 

■5s 



— 




(9b) 



To proceed further requires that the exact nature of be known. 

In the fully developed case^ ^ is a simple function of jT alone, 
and is easily treated. In the case of developing flows, such as the 
one under investigation, there is no precise, known functional descrip- 
tion of. 






B. DEVELOPING FLOW STABILITY DETERMINATIONS 

Schiller [2], and Langhaar [3] present workable models of developing 
velocity profiles in parametric form. However, after long and fruitless 
attempts, the author has concluded that there is no simplified mathematical 
solution suitable for representation of the developing mean flow at pre- 
sent. Accordingly, following Tatsumi [1], purely parallel flows which 
have the same velocity profile at a given axial distance as the developing 
flow shall be treated. For such velocity profiles, it was found that 
experimental data could be approximated to an excellent degree of 
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accuracy by a third degree polynomial obtained from a least-squares 
curve fitting procedure* To make the solution more general, it is 
assumed that 

from which 

:= ^ ( 11 ) 

is obtained* By obtaining the necessary terms from (11) and (6) and 
substituting them into (9), the following equations are obtained: 



where ^ “ -f ' ^ ^ (12b) 

and primes denote differentiation with respect to Y' . Treating (12a) 
first, it is clear that for large R, the terms on the right side of the 
equation will predominate, since even when is assumed small,ft(R is 
Still large* Yet, if the entire term is ignored, it is impossible 

to satisfy all of the^ boundary conditions of the problem. Order of 
magnitude considerations allow for the retention of the g** term, and 
thus satisfaction of all boundary conditions* Returning to the defini- 
tion of g, 






(13) 



Inserting (13) into (12b) then yields, 
where 

^ C<ii ( C-bo - b/ 

and 
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Before simplification, equation (14) is of the general form, 

y'V Pci-)'^'" i- Qir) /?W + To-) ^ = o 



a form which possesses a regular singular at the origin, r = 0* This 
fact guarantees the existence of at least one series solution of the 



form 






(15) 



which will converge for 0 r < r , where r is the distance from the 

o o 

origin to the next nearest singular point. As the disturbances intro- 
duced are confined to the vicinity of the pipe wall, r ^ 1, and again 
following Schlichting [6], the terms 




are neglected as they are all of the order 1, whereas the o^R term is 
assumed quite large. The resultant equation is 



Assuming a solution of the form (15), differentiating and inserting 



terms. 
















^5.0 J (17) 



^ L c>C 



R 












5 ^ 






o 
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Setting the coefficient of the highest power of r equal to zero results 



in the indicial equation. 






with indicies n = 0, 1,2,3 to guarantee that 0* 

Before proceeding further, specific attention is due to the boundary 
conditions of the problem* The no slip condition at the wall guarantees 



V* must be symmetric, hence zero on the axis, or f(0) = 0. Finally, 



= 0. Here, it should be noted that most previous works also include 
a condition that u*(0) = 0, or that f'(O) = 0. However, as noted by 
Leite [4], there is no guarantee that only those disturbances generated 
by the experimenter are being accepted by the flow, i«e., there is no a 
priori reason to expect u®(0) - 0* It is noteworthy that it is not neces 
sary to impose this artificial condition in this analysis. 

Taking the largest exponent of r first, it is found that the second 
highest ordered coefficient. 



efficient in turn, it is found that the recurrence relation yields: 



that f’(l) = 0, 




Since the flow is axially symmetric. 







(18) 



can only be identically satisfied of = 0. Continuing, for each co- 




(19) 




( 20 ) 
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a. = 



„ C <K a Q.0 









( 21 ) 



After substituting n = 3, then; 

a;, = 



^ 60 i(^(C-fcc) 



a<^ = 



?H0 






c »< R ( S fci+ 



With these coefficients, the final form of the solution becomes. 



■f ( r) “ (Xo i + ^ ^ C ^R. hi 






(22) 









The existence of additional linearly independent solutions is deter- 
mined by returning to the indicial equation. The value n = 0 will yield 
such a solution, but as the boundary condition f(0) = 0 would not be 
satisfied without a^ = 0, then this solution may be discarded. The same 
consideration applies for the value n = 1, i.e., the f®(l) = 0 boundary 

condition would not be satisfied without a = 0. In each of these cases, 

o 

if a^ were to be equal to zero, then it is evident that the series could 
be represented as 

^ /n + K+i 

> r 

T^o 

Therein, the values n = 0 and n = 1 would not then satisfy the resultant 
indicial equation. 

Finally, the case of n = 2, the last of the four index values is 
considered. Now, n = 2 identically satisfies the recurrance relation 
involving a^, i.e., the second highest ordered term. Thus, it is 
possible that this value would yield a linearly independent solution. 
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However, as this exponent differs from the highest by a positive integer 

(vis* 1), it is always possible to generate another linearly independent 

solution of the form^, 
oc> 

K-o 

This solution clearly fails the f(0) = 0 boundary condition, hence it is 
also discardedo Therefore, of the four linearly independent solutions 
possibly generated by the exponents from (18), only one remains that 
satisfies all the conditions of the problem, viz«, the solution represented 
in relation (22) o 



Applying first the wall boundary condition to this solution yields, 

-Tt-) = 1 - a ^ [(q.-w'-- J 

‘/c' 90 ?uo ^ ’ 

co(^tC,(c^-bo) ^ O 

Separating this equation into real and imaginary parts gives 



Cf2 ~ ^0 



/I - 



3 hz. — 

S'i/0 S'¥0 

84 0 P - ' yyo 



"(G-tol'-c’,'] 



= 0 (24a) 
= 0 (24b) 



Solving (24a) for C. results in 

1 



= -- n.¥^ (25) 



^Hildebrand, F. B., Operational Calculus for Applications , p. 137, 
Prentice “Hal 1 , 1962 o 
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Following the same procedure as above in applying the boundary condition 
f*(l) = 0 to the solution in relation (22) yields two additional equa- 
tions 






a 






hi hL _ ^ ^ ^ /'26a'» 

r VO /ao /CO ^ ^ 



^ . oC^RCc 

^ 8 ^ /:>sO 






Solving (26a) for gives 

(2'- ^ 32 h i ^ 



(27) 



It should be noted that C, and C have each been normalized with respect 

i r 

to U and that ^ has been normalized with respect to d* It is clear 
max 

that the normalized C value, the wave speed, takes on values 0 C < 1 
with the maximum value occurring in the region of maximum flow velocity. 
Inspection of (25) and (27) reveals that, within the developing region, 
must always be positive. Hence the developing region is stable. 
Further, in the limiting case, viz., in the fully developed flow regime, 
the coefficients of the polynomial representation of the velocity pro- 
file become; 



"o = 1. »3 ' -I- »2 = ^ ■ 0 



Analagous considerations demonstrate then that in the limiting case of 
fully developed flow, the solution represented by (22) also holds. ^ 



See Also reference 



[5]. 
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Data obtained from the curve fitting procedure is presented later 
to verify the applicability of this solution to the developing flow 
regime. The solution of (24b) would follow from Pekeris [5] under the 
assumption of purely parallel flows. Further pursuit of the exact solu- 
tion of (24b) however, depends on the exact determination of v as a 
function of r and z. 
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IIIo EXPERIMENTAL ANALYSIS 



A. TUNNEL FACILITIES AND EXPERIMENTAL PROCEDURE 

The experiments were conducted in a horizontal Lucite pipe 1.5 
inches I.D. and 45 feet (300 diameters) long. A centrigugal fan was 
used to pull air from the stagnation conditions in the plenum chamber 
through the system. The inlet chamber consisted of a rectangular box 
3 feet square in cross-section and 4 feet long which was fitted internal- 
ly with 6-24 mesh wire screens, followed by a 2 foot long settling cham- 
ber. Sizing of the chamber was fixed to limit the maximum interior 
velocity to 1.5 fps. The well rounded inlet horn was a quarter-section 
of a 9 inch by 18 inch ellipse molded of fiber glass and polished smooth. 
All joints were hand fitted to insure smoothness and uniformity of mating 
interior surfaces. 

The disturbance generator. Figure 1, was located 12 inches (8 dia- 
meters) from the pipe entrance, and was followed by measuring ports at 
16, 20, 32, 40, 48, 56, and 272 diameters downstream of the entrance. 
Uniformity of the internal diameter and concentricity of the pipe were 
found to be within .005 inch for all sections of the pipe. The indivi- 
dual pipe sections were joined by flanges, each pair hand fitted and 
joined to insure uniformity. The pipe sections joined to the distur- 
bance generator were fitted with rubber 0-rings to assist in damping any 
vibrations which might carry over from the driver to the pipe. The pipe 
was fastened to a horizontal metal track which was bolted to the con- 
crete deck. The vibrator was securely fastened to the deck after being 
isolated with one inch of corrugated rubber vibration mounting. The 
entire apparatus was longitudinally aligned with a high power transit 
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and leveled. Maximum deviation from true horizontal and longitudinal 
alignment is felt to be no more than .25 inch over the entire length. 

Measurements were taken with a Thermo -Systems constant temperature 
hot wire anemometer utilizing a boundary -layer type probe with a .095 inch 
long, .00015 inch diameter tungsten element. The probe could be located at 
any of the above mentioned stations, and accurately positioned to within 
.050 inch of the pipe wall. 

Prior to the commencement of data -taking, a careful survey was taken 
to determine the physical limitations of the apparatus. It was found 
that the maximum attainable Reynolds number was 9650; however, at this 
setting, the sensitivity of the flow regime was such that any disturbance 
within the laboratory was sufficient to cause transition of the flow. 

The maximum usable Reynolds number was found to be 8500. The survey fur- 
ther revealed that above a frequency of oscillation of 30 Hz., too much 
vibration was being carried over to the conduit. Hence data runs were 
limited to oscillations at 30 Hz. and less. 

The following procedure was used in taking data during the experi- 
mental phase: 

1. The mean velocity corresponding to the desired Reynolds 
number was measured and monitored at the station located in the fully 
developed flow region (272 diameters from the entrance) by regulating 
the amount of air pulled by the fan. 

2. The probe was positioned at the station nearest the genera- 
tor, and rotated until the position fixing the maximum reading was found. 

3. A velocity profile was first obtained; data taking runs 
then commenced with oscillation frequencies of 30 Hz., with % and ^ inch 
amplitudes, oscillations then at 20 Hz., with % and \ inch amplitudes, 
and finally at 10 Hz at % inch amplitude. 
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4. The probe was then moved to the next location down stream, 
and steps 2 and 3 repeated* 

5. Upon completion of recordings at all stations in the entry 
length region at the desired Reynolds number, step 1 followed by 2, 3 and 
4 were again repeated for another Reynolds number. 

The following procedures were used in recording data: 

1. The output of the hot wire was displayed on an RMS meter 
and on a dual beam oscilloscope, as was the signal input to the vibrator. 

2. As the probe was moved from its closest proximity to the 
wall towards the center of the tube in .025 inch increments, the hot 
wire bridge voltage and RMS voltage were recorded and the oscillogram 
either described or photographed. This step was repeated at each loca- 
tion for each data run conducted in (3) above. 

The data obtained in the manner just described was normalized 
through the use of appropriate values. Graphical presentations of this 
data are discussed in the following sections of this report. 

B. TUNNEL PERFORMANCE AND EXPERIMENTAL RESULTS 

1 . Developing Velocity Profiles 

Figures 2-5 show the velocity profiles obtained both at high 
and low Reynolds numbers at two points in the developing flow. They are 
compared with the analytical approximations of Langhaar [2], with a cor- 
rection factor of .8 applied to the development length. The correction 
factor was introduced because, similarly to Leite [3], it was found that 
the peak velocity obtained at a given axial location was always higher 
than that predicted by theory. It is additionally noted that the boundary 
layer itself does not develop as rapidly as predicted, although this 
factor is strongly dependent upon Reynolds number and thus may be a 
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characteristic of the apparatus. Figures 22 - 25 show the polynomial 
fitted by a least-squares technique, the results of which were used in 
the determination of the damping factor. The second curve in each of 
these figures is that of a fully developed profile offered for compari- 
son. Table 1 contains a listing of the coefficients of the polynomials 
used. 

The noise level in the developing flow is very much greater, by 
approximately an order of magnitude, than that previously reported in 
fully developed flow experiments. It is felt that due to the type of 
probe utilized and because of it*s orientation, some of the noise signal 
is indeed a measure of the radial component of velocity, v. In any case, 
the average turbulence level for all stations was .8 - 1.07o. As might 
be expected, the higher values occurred at close proximity to the wall, 
and the minimum value was found in the core region of the profile. 

The very nature of the hydrodynamic entry region, where most of the 
frictional losses are known to occur in pipe flow, suggest that one 
should expect a rather noisy signal. 

2. The Disturbance 

At the highest frequency and highest amplitude, the turbulence 
level jumped to a peak of 127 near the wall, with the peak at the center 
line of 8% at the farthest station located in the entry length. At the 
highest Reynolds number tested, the level at the station in the fully 
developed region reached a peak value of 1.07o compared with .87 without 
any disturbance. It was noted that at all Reynolds numbers and at fre- 
quencies and amplitudes accepted by the flow, that the radial propaga- 
tion of the disturbance was strongly dependent upon the velocity profile, 
i.e., upon the boundary layer thickness. Penetration into the core region 
was strongly amplitude dependent. 
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a. Reynolds Number Dependence 

One would expect that, for a fixed frequency and amplitude 
the level of turbulence at any downstream position along the pipe would 
increase with increasing Reynolds number* Figures 6-9 show plots of 
Turbulence Level vs* r/a at near and far stations at high and low 
Reynolds numbers and demonstrate this relation quite well. Figure 10 
shows the confining effect of the thinner boundary layer found at higher 
Reynolds numbers, i.e., the higher level disturbances are not allowed 
to diffuse towards the center into the core region. This effect de- 
creases dramatically with decreasing Reynolds number. Not only is diffus 
ion allowed, but it occurs sooner, as the lower mean velocity would 
indicate. Lastly it is noted that there were combinations of frequencies 
and amplitudes which produced extremely small fluctuations in the mean 
flow, just above the background noise level. Such phenomenon were 
observed at 10 Hz. at an amplitude of % inch and at 20 Hz. at % inch at 
all Reynolds numbers. It appears that these disturbances were simply 
damped within the boundary layer without producing the characteristic 
turbulent spike near the wall. 

b. Frequency Dependence 

Figures 11 - 14 show the very pronounced dependence of 
turbulence level upon the frequency of the disturbance generator. It is 
noted that as frequency is lowered, the speed at which the wave travels 
down the pipe should decrease. Hence the station at which a particular 
level of turbulence occurs for a given amplitude of disturbance at a 
given Reynolds number should be located further up stream as disturbance 
frequency is lowered. This is clearly shown in Figure 15. The idea of 
wave speed varying directly with the excitation frequency would seem to 
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suggest that, since the low frequency disturbance spends more time within 
the boundary layer, it interacts for a longer period with the strong 
viscous forces acting there. Hence it is damped more quickly and com- 
pletely. Again the observation that the velocity profile itself strongly 
effects the characteristics of the disturbance in the fluid is apparent, 
c. Amplitude Dependence 

Figures 16 - 19 show the increase in turbulence with 
amplitude at low and high Reynolds numbers# Figure 20 shows the strong 
dependence of diffusion upon amplitude, i.e., at higher amplitudes, for 
a given Reynolds number and frequency of disturbance, the higher the 
amplitude of disturbance, the greater the penetration of the disturbance 
into the core. Lastly, the fact that damping occurs at the frequencies 
and amplitudes tested suggests that there is an amplitude below which all 
traces of the disturbance would damp over a distance of a few diameters# 
Such a determination was not made here, but its existance is foreshadowed. 

3. The Damping Factor 

Following Leite [4], it is seen that from the definition of the 
u component of velocity and from (6a), 




If a suitable starting point, z , is chosen, this relation becomes, 

o 



or. 








- r 

iho') 



(29) 
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Thus , 



s Cp ^ 

c< ^ 



(31) 
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If values of ^ are plotted vs. the downstream displacement 

o 

(z - z )/d, the slopes of these curves will be, 
o 




Then, if either a value of wave speed, C^, or a value of wave number, CC 
can be determined, then can easily be computed. As a check, equations 
25 and 26 can then be solved to yield comparative values. Figure 26 

u* 

shows a plot of — vs. (z - data taken at Reynolds 

o 

number of 7600 at the maximum frequency and two amplitudes. As one 
would expect, the slopes increase as the center of the pipe is approached 
Further, at an r/a of .9, the larger disturbance is damped at a faster 
rate. 



Equations 24 were solved utilizing several experimentally obtained 
velocity profiles. Results are presented in Figure 27, which portrays 
the characteristic "hairpin” neutral stability curve. Plots are also 
presented for values of .05 and «1 which verify the stability of the 
flows tested. 
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IV. CONCLUSIONS 



(1) Flow in the laminar transition region of a pipe of circular 
cross-section is stable to small, rotationally syinmetric (torsional) 
disturbances. 

(2) There exists a strong dependence of amplification factor and 
wave propagation speed upon the shape of the developing velocity profile. 

(3) There is a dependence of both wave speed, C^, and wave number, 
c<. j upon radial position, though evaluation of this relationship will 

require an accurate description of v(r). 

(4) Within the range of frequencies and amplitudes tested, there 
existed a peak turbulence level in excess of 107. very close to the wall 
of the pipe. However, disturbances of this magnitude never reached the 
core of the flow and transition of complete turbulence never occurred. 

(5) Low frequency disturbances are completely damped within the 
developing boundary layer. 

(6) The experimental results suggest the existence of a **critical 
layer*' wherein the disturbance is sufficiently strengthened to cause 
transition to impend. Without this sufficient amplification, the dis- 
turbance diffuses into the central flow region, while continually being 
damped. 

(7) Small disturbance theory per se plays only a small role in the 
overall description of the transition of pipe flows from the laminar to 
the turbulent state. Better understanding of this transition phenomenon 
can only be achieved through an analysis which would describe the spatial 
growth of a finite disturbance rather than the temporal growth or decay 
of an infinitesimal disturbance. However, even the latter, if amplified 
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by the flow regime, quickly reaches an amplitude beyond which lineariza- 
tion of the equations of motion is no longer possible. Hence, there is 
at present no mathematical method suitable for a truely comprehensive 
study of transition. It appears then that further experimental studies 
are the methods whereby more information can be gathered on this subject. 
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V . COMMENTS AND RECOMMENDATIONS 



The experimental apparatus proved functional, easily maintained and 
operated, and very reliable ♦ It is recommended that further pursuit of 
this experiment be preceeded by the following alterations: 

(1) Construction of a new entrance section in accordance with 
the design of Cohen and Ritchie [8] to facilitate flow developing at very 
high Reynolds numbers* 

(2) Extension of the pipe to allow fully developed flow to be 
achieved at Reynolds number of order 15,000* 

(3) Integration of the air supply into a totally closed system 
to remove the possibility of interference from air currents and tempera- 
ture variations within the laboratory. 

(4) Modification of the disturbance generator itself for use 

of frequencies above 30 Hz. and amplitudes of %” from the present mechani- 
cal design to an electromagnetic design to eliminate vibration problems. 

After improvements are made to the apparatus, prior to the resump- 
tion of experiments in stability, it is strongly recommended that an 
experimental study of developing pipe flow itself be undertaken to deter- 
mine : 

(1) The developing velocity components u and v as functions of 
the pertinent geometric parameters r and z. 

(2) An analytical or numerical solution of the pertinent Orr- 
Sommerfeld equation. 

It is the author's belief that the study of the stability problem 
itself utilizing the small disturbance theory will ultimately lead to a 
better understanding of the problem of beginning of transition from 
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laminar to turbulent flow. Such a study, however, must utilize analyti- 
cal approaches more compatible with experimental procedure than those 
presented in the past. 
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TABLE 1 



Coefficients of the polynomial approximations to the velocity profiles at 

z/d ~ 16 and z/d = 56 



Reynolds 



z/d 


number 


b 

o 




^2 




16 


8500 


.702197 


-.752768 


2.75123 


-2.63370 


16 


5700 


.733547 


-.657861 


2.64586 


-2.7294 


56 


8500 


.823472 


-.378369 


1.89644 


-2.34190 


56 


5700 


.87101 


-.601644 


2.70048 


-3.02695 
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FIG 1 d/sturban/ce: (^enerrtor 
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Fisure 2: Developing Velocity Profile; a/d 16, R 5700 



39 




r 



Figure 3: Developing Velocity Profile; z/d 56, R 5700 
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Figure 4: Developing Velocity Profile; z/d 16, R 8500 
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Flgxire 5: Developing Velocity Profile; z/d 56, R 8500 
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Figure 6; Radial Diatrlbution of Disturbance Amplitude; 
z/d 16, R 5700 
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Figure 7: Radial Distribution of Dlsturbanoe Amplitude; 
z/d 48, R 5700 
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Figure 8 ; Radial Distribution of Disturbance Amplitude; 
z/d 16 , R 8500 
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Figure 9: Radial Distribution of Disturbance Amplitude; 
z/d 48, R 8500 
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Figure 10: Radial Distribution of Disturbance Amplitude* 
z/d 56, R 5700, R 8500 ’ 
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Figure 11: Radial Distribution of Disturbance Amplitude; 
z/d 16, R 5700 
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Flguro 12: Radial Dlatrlbutlon Of Dlsturbanoa Amplitude; 
z/d 48, R 5700 
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Figure 13: Radial Distribution of Dlaturbance Amplitude; 
z/d 16, R 8500 
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Figure 14: Radial Distribution of Disturbance Amplitude; 
z/d 48, R 8500 
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Plgxir# 15: Radial Distribution of Dlsturbanos Amplitude* 
z/d 16, R 5700 
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Plgxire 16: Radial Distribution of Disturbance Amplitude* 
z/d 32, R 5700 ^ ’ 
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Flsure 17: Radial Diatrlbutlon of Dlaturbanoa Amplitude: 
z/d 16, R 5700 
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PlKure 18: Radial Dlatrlbutlon of Disturbance Amplitude; 
z/d 48, R 5700 
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Flgur* 19: Radial Distribution of Disturbance Amplitude: 
z/d 32, R 8500 
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Figure 20: Redial Distribution of Dlsturbanoe Amplitude; 
z/d 48, R 8500 
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Figure 21; Radi«.l Distribution of Disturbance Amplitude; 
z/d 32, R 8500 
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Flgur6 22* Polynoniiftl Approxlmatrlon to tho Volocity Prof 11# } 
z/d 16, R 5700 
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Figure 23: Polynomial Approximation to the Velocity Profile; 
z/d 56, R 5700 
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Flgiire 24: Polynomial Approximation to the Velocity Profile; 
z/d 16, R 8500 
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Figure 25: Polynomial Approximation to the Velocity Profile 
z/d 56, R 8500 
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Figure 26: Decay of Disturbance Amplitude » R 7^00 
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